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Enhancing quantum transport in a photonic network using controllable decoherence 
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Transport phenomena on a quantum scale appear in a variety of systems, ranging from photo¬ 
synthetic complexes to engineered quantum devices. It has been predicted that the efficiency of 
quantum transport can be enhanced through dynamic interaction between the system and a noisy 
environment. We report the first experimental demonstration of such environment-assisted quan¬ 
tum transport, using an engineered network of laser-written waveguides, with relative energies and 
inter-waveguide couplings tailored to yield the desired Hamiltonian. Controllable decoherence is 
simulated via broadening the bandwidth of the input illumination, yielding a significant increase in 
transport efficiency relative to the narrowband case. We show integrated optics to be suitable for 
simulating specific target Hamiltonians as well as open quantum systems with controllable loss and 
decoherence. 


Recent research into photosynthetic antenna com¬ 
plexes has shown evidence of coherence in excitonic en¬ 
ergy transport HHa, despite the noisy cellular environ¬ 
ment in which such complexes are found. Indeed, envi¬ 
ronmental decoherence has been credited with increas¬ 
ing the efficiency of transport through these systems, 
an effect known as environment-assisted quantum trans¬ 
port (ENAQT) [5] or dephasing-assisted transport [3[Zl. 
While ENAQT has been the subject of many theoretical 
studies—whether in the photosynthetic context [EHle] or 
in other nanoscale transport systems —and de¬ 

spite its potential for improving transport in artificial 
quantum systems, it has so far never been directly ob¬ 
served. 

We use an integrated photonic simulator to demon¬ 
strate the first implementation of ENAQT. Our simula¬ 
tor was fabricated using femtosecond-laser direct writing, 
which allows waveguides to be drawn directly into glass 
using a focused pulsed laser. This permits the creation of 
three-dimensional waveguide arrays, as well as precision 
and repeatability in engineering interactions [2TH27]. We 
used control over the wavelength and bandwidth of the 
guided light to simulate effective decoherence, thereby 
enhancing transport efficiency. 

We consider a single excitation on a network of N cou¬ 
pled sites, governed by a tight-binding Hamiltonian [5] 

N N 

^=E Em \m){m\ + E ^mn {\m){n\ + \n){m \), (I) 

m=l n<m 

where \m) denotes the excitation being localised at site 
m, Em the energy of that site, and Vmn the coupling 
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between sites m and n. Although ENAQT can occur on 
an ordered lattice where all the energies Sm are equal [29] , 
transport enhancement was first explained in disordered 
systems, which we consider here. 

We are interested in the efficiency of transport from 
a particular initial site to a particular target site, where 
the excitation is trapped. In the case of a photosynthetic 
complex (figure [^), trapping describes the transfer of 
excitons to a reaction center, where they drive charge 
separation. It can be modelled as irreversible coupling 
of the target site to a sink at rate k, (figure [^). The 
efficiency is then the probability of finding the exciton in 
the sink after some particular time. 

ENAQT occurs when decoherence increases the trap¬ 
ping probability over the fully coherent case. Decoher¬ 
ence results from coupling of a quantum system to inac¬ 
cessible degrees of freedom. Eor example, in photosyn¬ 
thetic antenna complexes, the energies of chromophores 
are coupled to molecular vibrations; tracing out this envi¬ 
ronment results in decoherence in the excitonic subspace. 

In the absence of decoherence, energetic disorder tends 
to localise the wavepacket through processes such as de¬ 
structive interference or Anderson localisation [30H32], 
thus preventing it from reaching the target. Since these 
are coherent processes, they are diminished by decoher¬ 
ence, possibly resulting in enhancement of transport ef¬ 
ficiency to the target site. An alternative but ultimately 
equivalent view of ENAQT is that eigenstates of H are 
stationary, making it difficult to reach the target if the 
initial and target sites differ in energy. Incoherent pro¬ 
cesses, however, permit transitions between the eigen¬ 
states of H, yielding greater mobility. 

The first theoretical explorations of ENAQT focused 
on the case where the decoherence takes the form of site- 
independent, Markovian, pure dephasing I3IS]. Although 
the decoherence in our simulation is neither site-based 
nor Markovian, this is not an obstacle to ENAQT; [33113^ 
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FIG. 1. Environment-assisted quantum transport, a Photosynthetic antenna complexes are networks of chlorophylls that 
collect and transfer solar energy. A well-studied example is the Fenna-Matthews-Olson complex of green sulphur bacteria, here 
depicted as a network of 7 sites that transports excitation energy from initial site 1 to target site 3 (adapted with permission 
from ref. [281). Simulations have suggested that this transport may be enhanced by decoherence IMZl- b We simulate an 
instance of ENAQT on a lattice of 4 sites, with site 1 initially excited and site 3 the target. If the detuning of site 4 
equals C, one of the system eigenmodes has no occupancy at site 3 and cannot couple to the sink; by broadening the levels, 
decoherence breaks the condition A/3 = C, allowing all eigenmodes to couple to the sink and thus increasing transport efficiency, 
c Our simulator consists of 4 coupled waveguides arranged as shown (cross-section). The sink is modelled with a large array 
of closely coupled waveguides that transport light away from the main 4 waveguides. At the central wavelength Ao, waveguide 
4 has propagation constant /3 + A^d, while the others have propagation constant /3. d Theoretical expectation of transport 
enhancement for this system, as a function of simulation length z and decoherence strength 7 . The red bar indicates the region 
explored experimentally. 


transport efficiency can be enhanced as long as the de¬ 
coherence allows population transfer between otherwise¬ 
stationary eigenstates. 

We simulate ENAQT in an array of coupled single¬ 
mode optical waveguides obeying the equation USES] 

d ^ 

~ “I" 3 y C^mnO'hiz), (2) 

n^m 

where the light is propagating in the 2 ;-direction, is 

a creation operator for a photon in waveguide m at posi¬ 
tion and Pm and Cmn are respectively the propagation 
constants of the waveguides and the couplings between 
them. The former are determined by the waveguides’ re¬ 
fractive index profiles, while the latter also depend on the 
separations between them. Light propagation governed 
by this Schrodinger-like equation directly simulates evo¬ 
lution under iL, with Cmn replacing Vmn and the Pm 
replacing Cm- We can thus simulate different Hamiltoni¬ 
ans by controlling the number, position, and refractive 
indices of the waveguides. 

The intrinsic stability of laser-written waveguide ar¬ 
rays renders decoherence challenging to simulate. One 
approach is to stochastically modulate the index of re¬ 
fraction along each waveguide m Though every real¬ 
ization with a particular longitudinal index profile will 
be fully coherent, decoherence can be simulated by aver¬ 
aging over the recorded optical outputs from many real¬ 
izations in post-processing. This approach was recently 


used to simulate decoherence-enhanced navigation of a 
maze [38]. 

By contrast, we simulate decoherence by averaging 
over the results from a single array illuminated with many 
optical wavelengths. Although each individual wave¬ 
length propagates through the waveguide array coher¬ 
ently, decoherence can be achieved using broadband il¬ 
lumination and a single output intensity measurement 
which does not resolve wavelength. In other words, the 
wavelength degree of freedom is traced out yielding a 
partially mixed state. Similar approaches are well es¬ 
tablished in optics, where for instance thick birefringent 
quartz plates followed by wavelength-insensitive mea¬ 
surements have been used to decohere the polarization 
states of single photons [Ml HQ]. 

As an example of our approach to simulating deco¬ 
herence, we consider two uncoupled waveguides a and b 
that have a difference A/3 in their propagation constants 
at a wavelength Aq . As light propagates along the waveg¬ 
uides, it will accumulate a phase difference zAp between 
them. Broadband illumination can then be seen to cause 
effective decoherence on a length scale comparable to the 
illumination coherence length: Any initial coherence pab 
between the waveguides will decay: 

Pabiz) = Pab(O) 5^') > (3) 

where is the normalised first-order temporal 

correlation function of the light, which—being pro- 
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FIG. 2. Experimental setup and waveguide design, a Predicted dynamics of light in our waveguide array, as a function of z. 
The inset shows the predicted device output distribution for input at wavelength Aq. The sink array is sufficiently long that 
light reflecting from the far boundary fails to couple back into the system waveguides during the simulation, b A hbre-coupled, 
tuneable Tiisapphire laser in quasi-cw mode undergoes polarisation control (Pol) before it is imaged into the sample using a 
15 mm focal-length aspheric lens. The output is imaged via a 14x telescope onto a variable slit, which collectively measures 
the total intensity output from the system, bath, or all the waveguides using a large-area power-meter (PM). Alternatively, 
the output can be imaged onto a CCD camera for alignment and diagnostics, c CCD images of the output after 15 cm when 
illuminated at Ao = 792.5 nm and at 830 nm. As designed, the light in the system waveguides is evenly distributed at Ao apart 
from the target site, which is dark. In contrast, the target site is much brighter at 830 nm, indicating that more light will 
couple into the sink—a sign of ENAQT. 


portional to the Fourier transform of the spectrum— 
decays to zero faster for spectrally broader illumina¬ 
tion. The strength 7 of the decoherence can be quan¬ 
tified as the inverse of the optical coherence length, 

7 = (^^ iZo drj , and is usually proper- 

tional to the FWHM bandwidth A A. For a uniform dis¬ 
tribution centered at Aq, 


A/3AA 

27rAo 


(4) 


In the case of coupled waveguides, wavelength depen¬ 
dence affects couplings in addition to the propagation 
constants. The resulting decoherence will therefore not 
only have characteristics of pure dephasing (as in equa¬ 
tion (§), but will also include off-diagonal terms. How¬ 
ever, this kind of decoherence can also result in ENAQT. 
Indeed, the off-diagonal terms imply that our decoher¬ 
ence could, in principle, be used to simulate ENAQT in 
ordered systems [29], but the effect would be weaker be¬ 
cause of the absence of the pure-dephasing contribution. 

For our simulation we chose a network of four waveg¬ 
uide sites—arranged as in figure —because it is one of 
the smallest systems in which ENAQT is possible and 
because it can give significant enhancements even with 
relatively weak decoherence. Waveguide 1 was the in¬ 
put, waveguide 3 the target, and the sink consisted of a 
long linear array of tightly coupled waveguides. The cou¬ 


pling between waveguides in this sink was significantly 
higher than between the four main waveguides, so that 
any light entering the sink from waveguide 3 was largely 
transported away [Hjsa]- The sink need only be long 
enough to prevent light reflected from the far end from 
returning into the main simulation waveguides. 

We need to choose propagation constants and separa¬ 
tions among the four main waveguides in order to best 
approximate this Hamiltonian: 


fpco 0 \ 

C 13 C 0 

0 c ^ c 

\0 0 C /3 + A/3/ 


(5) 


where all parameters depend on the wavelength A of the 
input light. Due to their wide separations, couplings 
between non-neighbouring sites are negligible (< 5% of 
neighbouring-site couplings). Our simulator is designed 
so that, at our central simulation wavelength Aq, 


A/3(Ao) = C(Ao). (6) 

In this case, one of the eigenstates of H has no support 
on site 3, \ 2 pi) = (—1, —1, 0, l)/v^, while the remaining 
three eigenstates all have substantial support on site 3. 
Because cannot couple to the sink—at least at Aq— 
the maximum trapping efficiency at infinite time is 7 = 
l-(l|^i)=2/3. 
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Considering the wavelength dependence of H provides 
a different way to see the decohering effects of broad¬ 
band illumination. Due to the dependence of P and C 
on the wavelength, equation only holds at a particu¬ 
lar wavelength Aq. At other wavelengths, \p)i) will have 
some support on waveguide 3 and thus be able to couple 
to the sink, increasing the efficiency above 2/3, as shown 
in figure [^. Unlike in other examples of ENAQT [5], 
the efficiency increases monotonically with the strength 
of decoherence, meaning there is no optimal level of de¬ 
coherence in this model. 

Based on measurements of couplings and propagation 
constants in isolated pairs of waveguides (see Methods), 
we selected the following design parameters: Ap = C = 
1.0 cm“^, Ctrap = 1-5 cm“^, and Cgink = 1-75 cm“^ (see 
figure]^). Figure]^ shows numerical modelling of light 
propagation given these parameters. Although these 
were designed for a center wavelength of 800 nm, varia¬ 
tions in the implementation of the waveguide parameters 
resulted in equation being satisfied at Aq = 792.5 nm. 

We measured the efficiency using narrowband light 
(less than 1 nm bandwidth and always horizontally po¬ 
larised for consistency) from a tuneable Tiisapphire laser 
(Spectra-Physics Tsunami) in quasi-cw mode (figure [^). 
The output was imaged using a custom-built 14 x magni¬ 
fying telescope and the optical power was measured using 
a large-area power-meter after isolating either the system 
or sink waveguides using a variable slit. Examples of the 
output distribution are given in figure H, showing the 
significant difference between illumination at Aq and an 
off-center wavelength. Figure shows the measured ef¬ 
ficiency (fraction of light output in the sink modes) for 
wavelengths ranging from 745 to 835 nm. 

ENAQT—shown in figure —is the average enhance¬ 
ment in efficiency over the spectral band of interest, rel¬ 
ative to the efficiency at Aq, 

ENAQT = ~ 'n{>^o) ^ 

where the average (• • •) over A is taken over the top-hat 
distribution on (Aq — ^AA, Aq + ^AA). As predicted, it is 
an increasing function of the bandwidth, i.e., of the de¬ 
coherence strength 7 . The highest ENAQT observed was 
(7.6 ± 1 . 2 )% at a bandwidth of 95 nm. We thus demon¬ 
strated that coherent transport in a coupled, statically 
disordered system can be enhanced through decoherence. 

The theoretical prediction in figure contains no free 
parameters. It is a simulation of the dynamics under the 
simulated Hamiltonian i7, together with trapping from 
site 3 at the rate k discussed in Methods and pure de¬ 
phasing between waveguide 4 and the other three waveg¬ 
uides at rate 7 . The disagreement between theory and 
experiment is small considering the number of possible 
contributing factors. These include the off-diagonal de¬ 
coherence when the waveguides are coupled, the fact that 
the trapping is not perfectly exponential, errors in the 


measurements of the coupling constants, optical losses, 
and error in satisfying equation As an example, the 
shaded band in figure represents the error that would 
arise if AP{Xq) deviated from C(Ao) by up to 10%. 

In our experiment, the magnitude of ENAQT was lim¬ 
ited by the maximum achievable decoherence, which at 
7 = 0.02 cm“^ was small compared to the inverse of 
the propagation length. We were limited by two com¬ 
ponents of equation Q: the tunability of our laser lim¬ 
ited A A, while Ap was limited (via equation by the 
need to keep C small enough to stay in the tight-binding 
approximation. These are not fundamental limitations, 
and future improvements that increase 7 would result in 
significantly larger transport enhancement; as shown in 
figure [^, the transport efficiency of this model can get 
arbitrarily close to 1 for sufficiently long propagation dis¬ 
tances and decoherence strengths. Stronger decoherence 
would also allow ENAQT to be observed in networks that 
are less sensitive to decoherence than our model. 

Our results demonstrate that integrated photonics are 
well-suited for simulating open quantum systems, capa¬ 
ble of implementing a disordered target Hamiltonian with 
controllable loss. The technique of using broadband exci¬ 
tation to introduce tuneable levels of decoherence ensures 
that photonics can simulate open quantum systems; pre¬ 
viously, low intrinsic noise in photonic devices rendered 
decoherence difficult to realize in integrated optics, par¬ 
ticularly without averaging over results from many dif¬ 
ferent device realizations. Our approach will not only 
improve the controllability of photonic quantum simu¬ 
lation, but also aid in the experimental optimisation of 
transport in other engineered quantum systems. 

METHODS 

Our experimental setup is shown in figure The 
waveguides were fabricated in high-purity fused sil¬ 
ica (Corning 7980) using a laser direct-write technique 
whereby Tiisapphire laser pulses are tightly focused into 
the sample, which is then translated in three dimensions 
to yield continuous regions of positive refractive index 
change which act as waveguides [22]. To obtain single¬ 
mode waveguides with the desired propagation and cou¬ 
pling characteristics, laser pulses of duration 150 fs, en¬ 
ergy 400 pJ, central wavelength 800 nm, and repetition 
rate 100 kHz were focused 400 pm below the surface using 
a 40 X microscope objective and translated at 75 mm/min 
(for waveguides other than 4). 

Permanent index changes result in the focus, yield¬ 
ing elliptical, vertically oriented, single-mode waveguides 
with modes approximately 17x19 pm in size at 800 nm. 
Waveguide 3 and the sink waveguides are in a plane par¬ 
allel to the surface, while angles of 120 ° between this 
plane and the other system waveguides minimise non¬ 
neighbour coupling. Because couplings between elliptical 
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FIG. 3. Magnitude of observed ENAQT. a Transport efficiency—the portion of light that makes it to the sink—as a function of 
wavelength. The minimum efficiency at Ao = 792.5 nm is 77 = 0.636 ± 0 . 002 , slightly less than the theoretical inhnite-time limit 
of 2/3. The error bars are standard deviations caused by imperfect repeatability in coupling light into the sample and laser 
power fluctuations, b ENAQT—the relative increase in the efficiency over that at Ao—as a function of the optical bandwidth 
(top horizontal axis) and corresponding decoherence strength 7 (bottom horizontal axis). The red points are obtained by 
averaging the measured efficiencies over a uniform broadband spectrum with width A A and centered at Aq. The blue line 
represents the theoretical ENAQT, calculated based on the model in hgure and containing no free parameters (it is the 
cross-section of hgurealong the red line segment). The shaded region represents possible ENAQT if A/7(Ao) deviates from 
C(Ao) by up to 10%. 


waveguides are dependent on angular orientation [21] , we 
determined the couplings by writing pairs of waveguides 
oriented at the specihed angles. Couplings as a function 
of separation were determined by writing the pairs at 
different separations and measuring the output intensi¬ 
ties after a known propagation length when only one is 
optically excited. 

The hnal waveguide separations were chosen to ensure 
that the tight-binding approximation is maintained, that 
C can be matched by A/3, that Cgink > Ctrap > C', 
and that the number of sink modes remains manageable. 
The system-sink population transfer can be modelled as 
a constant effective rate [41] k = 2x^(1 — if 

X = Ctrap/Csink ^1- In our case, although x = 0.86 
meant that the decay was not perfectly exponential, it 
was nevertheless effectively irreversible, considering that 
sufficiently many sink modes were present to prevent re¬ 
flection from the far end and back into the main waveg¬ 
uides. 

For waveguide 4, the writing translation speed was de¬ 
creased to increase the propagation constant. In a pair 
of waveguides with mismatched propagation constants, 
the maximum power transfer depends only on the ratio 
Ap/C. We used this fact to determine the translation 
speed decrease necessary to set A/3 = C at 800 nm. 
Due to the measurement uncertainty (about 8%) and 
day-to-day variability in the waveguide writing process, 
several arrays were written with slightly different writ¬ 
ing speeds for waveguide 4 to ensure that one set would 
yield A/5 = (7 in our desired wavelength range. In the 


best sample, the writing speed was decreased by 12%. 
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